Relaxation of the Bose-condensate oscillations in the mesoscopic system at T=0 
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The general system is given of nonlinear equations describing dissipationless evolution of the 
oscillating Bose-condensate. The relaxation of transverse oscillations of the condensate in a trap of 
the cylindric symmetry is considered. The evolution occurs due to parametric resonance coupling 
the transverse mode with the longitudinal ones. The nonlinear rescattering in the subsystem of 
discrete longitudinal modes results in suppression of the return of energy, yielding dissipationless 
nonmonotonic relaxation of transverse oscillations in the condensate. 
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The discovery of Bose-condensation in ultracold gases 
has yielded a unique possibility for the study of evolution 
of the coherent properties in the macroscopic system iso- 
lated from environment. One of most interesting aspects 
here is connected with clarifying and analyzing the nature 
of the damping of coherent oscillations in the condensate. 
So far, the theoretical and experimental investigation of 
the problem has been connected, in fact, with considera- 
tion of the damping at finite temperatures (see detailed 
bibliography in IJ). In this case the ensemble of normal 
excitations plays effectively a role of the interior thermo- 
stat. The interaction with thermal excitations results in 
damping the condensate oscillations with the transfer of 
energy into the subsystem of normal excitations. 

A question about the internal mechanism of damping 
at T = represents a special interest in this problem. 
In the work of authors 2] the existence of such mecha- 
nism is shown on the example dealing with the damping 
of the radial coherent oscillations in the condensate with 
the cylindric symmetry in the transverse parabolic poten- 
tial. It is found that the damping arises from a special 
parametric resonance resulting in the transfer of energy 
into the subsystem of longitudinal modes. The paramet- 
ric resonance is due to oscillations of the sound velocity 
as a result of transverse oscillations of the condensate. It 
should be emphasized that the results obtained in can 
be used only for describing the initial period of damping. 

To describe the damping at the large times in the iso- 
lated system of finite sizes with the discrete system of 
energy levels, it is necessary to know additionally the 
temporal evolution in the subsystem of transverse modes 
and involve variation of the parametric resonance with 
reducing the amplitude of transverse oscillations. This 
aspect is a significant difference from the conventional 
picture of parametric resonance. The choice of geometry 
of the system is determined by the fact that, for an ar- 
bitrary variation of the frequency in a two-dimensional 
parabolic potential of the circular symmetry and arbi- 
trary magnitudes of parameters, as is strictly found in 
(also PI), two-dimensional oscillations do not damp. 

Recently, Paris group JJ has published the results of 
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the experimental study on the damping of transverse 
modes in the condensate as a whole (breathing mode, 
BM) in the elongated trap with the azimuthal symme- 
try for extremely low temperature of 40nK. At the small 
magnitude of the BM amplitude the authors observe the 
record slow damping. For larger magnitude of the initial 
BM amplitude, the picture changes drastically. After the 
fall within a limited time interval the reverse transfer of 
energy with the growth of the BM amplitude starts. The 
normal behavior of damping recovers only a noticeable 
time interval later. 

The theoretical analysis with explanation of the anor- 
mal picture is given in l5l|. It is found that the involve- 
ment of nonlinear coupling between longitudinal and 
transverse modes is essential. However, the relaxation 
of longitudinal excitations was considered phenomeno- 
logically in essence. 

In fact, for the discrete energy spectrum inherent in 
the closed mesoscopic system, the irreversible damping 
does not appear. It can be spoken only about the non- 
stationary redistribution of energy in the subsystem of 
longitudinal modes. 

In the present work we have obtained a full set of 
nonlinear equations describing the temporal evolution of 
the oscillating condensate in the lack of the irreversible 
channels of dissipation. Our starting point is the Gross- 
Pitaevskii equation. The condensate of the cylindrical 
symmetry is considered in a transverse parabolic poten- 
tial with longitudinal size 2L ^ R where R is the static 
radius of the condensate. The analysis of the solution 
demonstrates an essential dependence of the damping of 
radial oscillations upon the character of evolution of lon- 
gitudinal modes. Redistribution of the energy transferred 
into this subsystem results in the chaotic nonstationary 
picture of filling the discrete levels and simultaneously 
in the significant reduction for the amplitudes of active 
longitudinal modes excited directly as a result of para- 
metric resonance. This reduction is equivalent effectively 
to relaxation but occurs in the lack of real dissipation. 

We will consider the dynamics of Bose-condensate in a 
rarified gas at T = within the framework of the nonlin- 
ear Gross-Pitaevskii equation for condensate wavefunc- 
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+ [/o**(r,i)*'(^i)- (1) 

Here Uq = ATrti^a/m is the vertex of the local interaction 
between particles and a is the scattering length. The sys- 
tem is assumed to be closed and in the course of evolution 
the number of particles N and total energy E conserve. 
We restrict ourselves by the case of repulsion between 
particles (a > 0). Let us represent the wavefunction in 
the form 



The imaginary part of Eq. yields directly the equa- 
tion of continuity for the condensate 



dni {r, t) 
dt 



(2) 



v{f,t) = — V>I>(f,t). 
m 

In Q we introduced the fraction of density ni (r, t) de- 
pending explicitly on the time and represented the total 
density as n (r, t) = uq (r) + ni (r, t). 

The real part of Eq. |^ yields an equation for the 
phase 
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Below we consider a gas of the sufficiently high density 
so that the correlation length ^ = 1 / ^/mnoUo is small 
compared with all sizes of the system. It proves to be 
that only the long wave condensate excitations with the 
wavelengths larger than ^ are involved into the processes 
considered. In these conditions corresponding to the 
Thomas-Fermi approximation the last term in ^ can be 
neglected, e.g., For the stationary case when v ~ 0, 
the phase of the condensate wavefunction has the familiar 
value $0 = —fJ-ot/fi where no is the chemical potential. 
Then from equation Q one has 



Within the linear approximation the system of equa- 
tions ^ and ||SJ) reduces to the equation of harmonic 
oscillations 



Here 



Join = -: : '^00 = no[0), Cp = _ = — . (6) 



m m 



Let us introduce the whole orthonormalized system of 
eigenfuctions of the harmonic problem {Xsi'r)} being so- 
lution of the equation 



^IXsir) + c^V (/o(f)Vxs(r)) = 0. 



(7) 



Let us expand ni and $ in the whole system of eigen- 
functions Xs(f) 



ni 



{t,r)=Y,r.sit)Xsir), $ = 5]a,(i)x,(f) (8) 



Inserting these expansions into Eqs. Q and ©, one 
finds 
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Here (...) — J (Pr{...). For the right-hand side of the first 
equation, we employ the transformation 

(x:v(xs, (VxsJ)) = - ixs, (Vx:) (Vx.J) 

Let us introduce new variables 

For the modes described by real eigenfunctions Xs, coef- 
ficients Us, bs are real and y-s = Us- 

For these variables, Eqs. Q can be represented as 



no(f) 



(A'o - V{r)) 



(4) 



Keeping notation <I> for the phase connected with the 
condensate dynamics, we find 



-h— ^ — (V$f 

dt 2m^ ' 



UqUi 



(5) 



A set of nonlinear equations (0) and ^ describes not 
only excitations of quantum fluid but also their interac- 
tion. We will strictly take this interaction into account, 
not involving the perturbation theory. Note that the ex- 
ternal potential enters equations ^ and |SJl only implic- 
itly via the static distribution of density 
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Here 



^^Kss^sAvsi +y*-si)iys2 +y-s2) (lo) 
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Mss.s2 = 2k (S&j^'^ ii^xDXs. (VxsJ) , ^^^^ 
Kss.s2 = ^ (gfc)'^' ixt (Vx.J (Vx.J) . 

Let at the initial time moment a transverse condensate 
oscillation alone be excited as a whole with conservation 
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of the cylindrical symmetry ("breathing mode" and here- 
after index _L). A set of equations (|10|l can be represented 
in this case as 



k 

-i^K^k-kivk + ylk){y-k + yl), 

k 

^Vk + iVLkyk = ~'2Mkk±ylkiyi- + yl) 
+K±k-kiy± - y*i){yk + ylk) + ^k- 



(12) 



In the second equation we use equalities Mk±k = 
K±_k-k, Kkk± — Mkk±- The last term on the right- 
hand side determines the evolution in the longitudinal 
subsystem of excitations 

h = ^Mkk^k2{yk^ -y-ki){yk2 +ylfej 

- \^Kkk,k2{yki + y-k^)iyk2 +y-k2)- (13) 

fclfe2 

The eigenfunction of the breathing mode is found from 
the solution of equation Q involving that within the 
Thomas- Fermi approximation fo — ^ — r^/R^ where 
= 2cq/w^, uj± being the frequency of a parabolic 
trap. 

X± = /|(1-|J), n^ = 2Lu^, (14) 

Quantity k determines the value of wave vector for longi- 
tudinal excitations, running discrete values due to finite 
sizes of the system in the z direction. Following work ^ , 
it is easily to find the eigenfunctions and eigenvalues for 
the longwave longitudinal modes 



Xk — y y 

flk — ck{l 



(kR) 



2r^ 



{kRY 



16 R 
) C = 



(15) 



Using H14|l and (|15|l . one can straightforwardly calculate 
matrix elements (lllfl entering in (|12|) 



M^k-k 



Ml 



kkl^ 



Ve ( Ofc 
''4 V <^ , 



K.k-k = -v'-^{^) , 
and in (O 
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(16) 



Mkkik2 

Kkkik2 = 



sign(fcfc2)(5, 



k,ki+k2 I 



£1 £1 \ 1/2 

" j sign(fcifc2)<5fe,fci+fc2 



(17) 



Equations H12|l and (|13|l with matrix elements H16|l and 
(|17|l describe fully the relaxation of the transverse con- 
densate oscillations in the mesoscopic system in the lack 
of dissipative channels for evolution. 



We assume that the initial amplitude of breathing 
mode is relatively small, i.e., \5R/R\q <C 1. This re- 
stricts a scale of the nonlinear interaction between modes. 
From the other side the energy interval in which the dis- 
crete longitudinal modes are connected effectively with 
the breathing mode proves to be small compared with 
uj±. Under conditions we can simplify a set of equations, 
restricting ourselves with the quasi resonance approxi- 
mation. Within the framework of the approximation in 
equations (|12|) and H13I) we can retain only the terms for 
which the following inequality is valid 



Ail = liis ± ± 



(18) 



Let us rewrite equations 1)12(1 and H13|) within this ap- 
proximation, substituting ratios j/s — ys{t)/ \y±{0)\ for 
amplitudes ys ■ The initial amplitude of breathing mode 
y±{0) can be found from comparing the vibrational 
energy of the condensate at the initial time moment 
Evib = ^fioN \SR/R\l with the energy of transverse mode 
2hni_ |y±(0)|^ Hence 

\y±m 
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SR 






R 



After involvement of these notations equations 112|l and 
((T^ go over into 



d_ 

di 



y± + iil±y± 



k>Q 



A dash in a sum of the first equation means summation 
over "active modes" alone, interacting directly with the 
BM according to restriction H18|l 
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—yk + lilkyk 
ot 



2a 



'E 

fcl>0 



ay-kyi- 

ilkilki^k-ki 



fc>fci>0 
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ykiyk-ki 



Qk^ki^k+ki 
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ylVk+k,, (fc>0) (19) 



For longitudinal modes which do not interact directly 
with the transverse mode, one should omit the first term 
on the right-hand side of Eq. (|19|l . In these equations 



R lo ■ 



a ^ \{2M^k-k - K^k-k)\\y±m 
a' « 2a. 

(20) 

Let us take that one has yfe(O) = y_A;(0) for the trivially 
degenerated states at the initial time moment. It follows 
from a set of equations (|19|l that this equality holds in the 
course of evolution. To analyze, it is convenient to select 
the fast phase from complex variables ys, representing 
them as 

y,=bs (t)e-("=*+^»(*)), 
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where bs (t) and (fs (t) are the real quantities varying rel- 
atively slow in time at a/Lu± <C 1. Let us introduce 
notations 

7fefeife2 = i^k - ^ki - ^k2)t + iVk - Vki ~ Vk2) (21) 

Separating real and imaginary parts H19|l . one finds 

d 



k>0 



feC0S7_Lfefc, -g^ip± 



k>0 



^2_, ^sm7_Lfcfc, 



(22) 



1,0 




0,8 


w — 




5=1,2 


0,6 
^ 0,4 




\ 5=0,7 


5=0,95 \, 


0,2 








0,0 


^=0.5 



= cos7j_fefe + 4, 

4 = a 2^ I Lj bk^bk-ki COS jkkik~ki 

fc>fci>0 

-2q; I ^1 Ofc+fciOfci C0S7fe+fcifcifc, 

fei>0 

(23) 



/I 



fc>fci>0 



-2a' 



fci>0 



I ) 

I — — J 



sin7^fc/c + I'k, 



Sm^kkik-ki 



■smjk+kikik- 



(24) 

The analysis of the equations obtained allows us to 
find a scenario of relaxation of the transverse condensate 
oscillations. At the initial time moment the longitudinal 
modes are not excited and relation bk (0) ^ 1 is valid for 
them, while according to definition b± (0) = 1. Thus at 
the initial time period the right-hand sides of Eq. H22() 
and terms Ik, I'k quadratic in bk play no role and evolu- 
tion, in essence, is governed by first terms on the right- 
hand sides of equations H23|l and (|24|) with constant b±. 
The joint solution demonstrates an exponential growth 
bk(t) provided 



\An\ 



nk 



< a . 



(25) 



In particular, for |Af2| = 0, one has "f±kk — and equa- 
tion (|^ yields bk{t) = bk (O)expai. These results, in- 
volving requirement of the finiteness for initial amplitude 
bk (0), are a typical manifestation of parametric reso- 
nance (cp. 01 )■ AH longitudinal modes within the en- 
ergy interval of about a {active modes) experience an 
exponential growth. For the length of cylinder 2L, the 
spacing between levels equals Sfl = nc/L. If Sil < 2a, 
at least, one mode lies within this interval. The growth 
of bk{t) results in reduction of b±{t) after some delay 
and transfer of energy to the other parallel modes at the 
same time. For L ^ R, there is a large number of such 
modes but a discrete character of energy spectrum results 
in the total prohibition of irreversible processes. How- 
ever, irregular character of the transition amplitudes like 



FIG. 1: The behavior for the population of breathing mode 
as a function of the time at fixed value ko = IOtt/L for ^ = 
0.5, 0.7, 0.95, and 1.2. 




5=1,2 



FIG. 2; The behavior for the population of active mode ko 
as a function of the time at the same values ^. The plots are 
shifted with respect to each other by 0.1 along the vertical. 



cos7fcfejfe2 ill and in nonhnear 

terms Ik, I'k leads to chaotic evolution in the system of 
longitudinal modes. One should think that in this case 
the return to the active longitudinal levels should sig- 
nificantly be suppressed and a noticeable fraction of the 
transferred energy should remain in inactive modes. 

The direct numerical simulation of system H22|l - H24|) 
displays the picture described. Let only single level 
ko lie within the energy interval of about a. We ne- 
glect weak dispersion of longitudinal modes l(T^ . Then 
"""jj""" ^ kk^ expression ^ for longitudi- 

nal modes one can omit the first term. Let us introduce 
dimensionless time t = at and take into account that 
a'/a = const (see, (j2I3))- Then the evolution described 
depends, practically, on ratio f — Afl/a alone and, to 
weak extent, on kg for koL/n ^ 1. (Note that this is 
valid not only in the quasi resonance approximation but 
also for the general system of equations 1)12(1 .^ In Fig.l 
the dependence b\{T) = {SR{t)/SR{0)) is given at vari- 
ous magnitudes of parameter ^ for fixed value kg (in units 
tt/L). Here we involved not only the discrete levels ly- 
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ing below the active level but also the longitudinal levels 
above fco which start to be occupied at the next stage af- 
ter decay of state ifep. The law of energy conservation, 
which system (I22ll - (|24|) satisfies, results in a weak popu- 
lation of the upper levels in the course of evolution. In 
calculations we restricted ourselves by the same number 
of levels above and below ^Iko- For value ^ = 1, 2, the 
transverse oscillations does not decay at all. This is obvi- 
ously seen from Fig. 2 in which the behavior of population 
b^^ of level fco on r is plotted for the same values of pa- 
rameter ^. For given value ^, magnitude b^^ (r) remains 
close to zero for all r. The statement that the damp- 
ing of oscillations at T = should be absent at all with 
violation of the conditions for appearance of paramet- 
ric resonance (|25ll finds thus a direct confirmation. The 
curves corresponding to ^ = 0.7 and 0.5, on the contrary, 
demonstrate an origination of damping with the nontriv- 
ial character (see. Figs. 1 and 2). Reduction b\{T) and 
growth 6^^^ (r) are followed by the reverse transfer of en- 
ergy to the BM and the decrease of the population of 
the active mode. Only later, when the redistribution of 
energy over the other longitudinal modes becomes no- 
ticeable, there arises a traditional monotonic behavior of 
relaxation. The energy, as direct calculations show, is 
spread over all longitudinal modes, their population ex- 
periencing chaotic evolution. At this stage, in spite of the 
lack of irreversible processes, the character of interference 
due to dispersion of dynamic phases prevents from any 
noticeable growth of amplitude 6fc(, (r) and, thus, from 
reverse transfer of energy to the breathing mode. Then 
the evolution looks like relaxation of the transverse con- 



densate oscillation with concentration of larger fraction 
of energy in the subsystem of longitudinal modes. 

For the first time, nonmonotonic character of the 
damping of transverse breathing mode, analogous to that 
shown in Fig.l, is observed experimentally in Large 
magnitude \6R/R\q provides condition ^ < 1. It should 
be noted that the theoretical results obtained are univer- 
sal to essential extent. Thus the displayed picture of the 
effective relaxation with its anomalously nonmonotonic 
behavior has a general character at T = 0. 

However, quantitative comparison with the experimen- 
tal results U is difficult, in the first turn, due to the 
finiteness of temperature. Though T < /xq in the exper- 
iment, but T > hu!± and the longitudinal levels prove 
to be temperature-occupied at the initial time moment. 
(Note that with the introduction of random dispersion of 
phases ipk at the initial time moment the qualitative pic- 
ture of evolution holds but quantitative picture changes. 
The slow damping close to monotonic and observed in 
at sufficiently smaller value \6R/R\fj has analog in an 
isolated system at T = only for close to unity (see, 
Fig.l and ^ = 0.95). However, one cannot exclude that 
such damping is associated with the dissipation processes 
due to external factors. For sufficiently large magnitude 
of ratio \SR/R\q, these processes cannot play a role at 
most interesting stages of nonmonotonic relaxation. 
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